Robustness of channel-adapted quantum error correction 
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A quantum channel models the interaction between the system we are interested in and its environ- 
ment. Such a model can capture the main features of the interaction but because of the complexity 
of the environment we can not assume that it is fully accurate. We study the robustness of quantum 
error correction operations against completely unexpected and subsequently undetermined type of 
channel uncertainties. We find that a channel-adapted optimal error correction operation does not 
only give the best possible channel fidelity but it is more robust against channel alterations than 
any other error correction operation. Our results are valid for Pauli channels and stabilizer codes 
but based on some numerical results we believe that very similar conclusions can be drawn also in 
the general case. 

PACS numbers: 03.67.Pp 



I. INTRODUCTION 

One of the central problems of the quantum informa- 
tion processing systems is to find suitable error correcting 
procedures. Therefore Shor's early result which was 
followed by the development of the standard theory of 
quantum error correction (QEC for short) 0,0,131, was 
a very important step, giving the possibility for quan- 
tum information theory to become a potentially appli- 
cable science from an only theoretically interesting field. 
These standard quantum error correcting methods tried 
to adapt the techniques of classical information theory 
for the case of quantum information. However, it was 
later recognized that looking for the best error correcting 
method can be formulated as an optimization problem 
[S B > and recently some authors approached the task 
of error correction from this new viewpoint [1, H, Oj], EU , 
and found that better solutions can be achieved than by 
the standard QEC. But there are relatively few studies 
about the comparison of these approaches in the aspect of 
resistance against the uncertainty of the channel model. 
In this paper we take a step in this direction. 

In QEC the Hilbert space of the system of interest is 
embedded into a Hilbert space of a larger system. This is 
the encoding. In the most simple case this means that a 
logical qubit is encoded into the strongly entangled states 
of a few physical qubits. After encoding the system goes 
through a noisy quantum channel then we apply a QEC 
operation to eliminate the effect of noise. The method 
of standard QEC is based on the assumption, that if the 
noise is weak, then those errors will be the most probable, 
which only corrupt the state of very few, perhaps only one 
qubit, therefore the standard error correction procedure 
is designed to correct perfectly these types of errors. This 
assumption is true under very general circumstances, and 
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it is independent from the actual character of the channel. 
We only need to necessitate that the noise is acting on the 
qubits independently, but nothing else about the specific 
form of the channel. 

However, the channel models the interactions of the 
physical system with its environment, about which we 
can make plausible assumptions based on knowledge 
about the specific system. If we have information on the 
specific effect of the channel, then based on this we can 
design noise-specific encoding and error correction. The 
statement of error correction in this context, as an op- 
timization problem on ly g ained significant attention in 
the last few years d, HI, OHII ■ They showed for example 
that by having knowledge about the exact effect of the 
channel, the error can be corrected with better efficiency, 
than by conventional QEC. A smaller system is enough 
for an equally effective encoding, and/or the code offers 
better protection. 

However it is also undoubted that in reality, because 
of the complexity of the environment, the interactions 
between the system and its environment can not be taken 
into account with perfect accuracy. The quantum channel 
will always remain just a model, which though describes 
reality more or less precisely, but we can not assume that 
it is fully accurate. So the question rises: if the noise 
differs in some degree from the chosen model, then how 
effective does our model-specific error correction method 
remain, i.e. how robust it is? 

In this work, our aim is the study of this question, that 
how robust are the conventional and optimization based 
channel specific error correction procedures against the 
uncertainties of the channel. At first sight we could think 
that the conventional QEC — which makes very slight as- 
sumptions about the properties of the channel — is more 
insensitive to the uncertainties, than an error correction 
optimized for a specific channel. On the other hand, we 
can also think that if a solution is optimal, then in gen- 
eral in case of the very small alteration of the conditions 
it decays only in second order, so it is robust in some 
measure. These questions apparently have greater signif- 
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icance when the channel is general, but because of the 
simplicity of the handling of Pauli channels, here we con- 
centrate mainly on this specific case. 

The paper is organized as follows. In section [II] we de- 
fine the measure of robustness. In section UIT1 we demon- 
strate the procedure of searching the optimal correction 
operation for Pauli channels. In section lTVl we present the 
general model of robustness for Pauli channels. Finally, in 
section [V] we make conclusions related to Pauli channels 
and based on numerical results we discuss the possible 
generalizations of these results to non-Pauli channels. 



II. ROBUSTNESS 

A quantum channel is defined to be a completely pos- 
itive and trace preserving (CPTP) map [12J. This map 
can be represented by a set of operator elements, {E,}, 
called the Kraus representation, which gives the channel 
output as 



(1) 



and the operator elements must satisfy the completeness 
relation 



i 



(2) 



in order to represent a trace preserving map. We have to 
mention that the set of operator elements is not unique. 

First of all we have to characterize channel perfor- 
mance. The role of the entangled states in storing and 
transferring quantum information was explored long ago, 
and it is well established that the entanglement fidelity 
properly measures how successfully the channel preserves 
the entanglement, i.e. the quantum information [13j . But 
the entanglement fidelity depends not only on the chan- 
nel £ , but also on the input state p. To avoid this, in this 
paper we always use the maximally entangled state as an 
input. The argument behind this is that if the channel can 
preserve the maximally entangled state, then it can pre- 
serve other entangled states also. This special case of the 
entanglement fidelity is sometimes called channel fidelity 
@, as it depends only on the properties of the channel. 
This can be given in terms of {E^} Kraus operators as 



F ch (f) = Ei Tr (^)i 2 = E 
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where TL is the Hilbert space of the system. 

How effectively does a QEC operation TZ correct the er- 
rors which occur while the system goes through channel 
£ ? As usual, in this paper we measure it with the perfor- 
mance of the composite channel lZo£, that is with the fi- 
delity F c \ l (lZo£ ). The aim of the optimization based chan- 
nel specific error correction is thus to determine the error 
correction channel operation TZ*, for which F^IZ* o £) 



is maximal, so the result certainly depends on the chan- 
nel. The standard QEC also makes assumption about the 
channel model, but much less, as we have already men- 
tioned in the introduction. 

By fixing an encoding and a channel, we can compare 
the two methods. As a first step, we can compare the effi- 
ciency of the different error corrections, measuring them 
in terms of channel fidelity, certainly getting the opti- 
mization method always as the best. This however may 
not be striking, as if we know the nature of the error 
source, we can act against it more effectively. 

As a second step, we can compare the robustness of 
error correction operations. We can approach robustness 
from different aspects depending on how we model chan- 
nel uncertainty. One possibility is to take the uncertainty 
into account in the formulation of the optimization prob- 
lem, i.e. define the channel as an integral weighted by 
a probability distribution, then maximize the average 
entanglement fidelity. This method is shortly discussed 
by Fletcher in [Toj ] . A similar problem statement can be 
found in the work [TTj . which tries to solve the problem 
with a reformulated, "indirect" version of the optimiza- 
tion problem. Of course, to use these approaches we must 
first consider by what assumptions the distribution will 
be physically realistic. 

The drawback of such approaches is the exactly de- 
fined, "well known" manner of uncertainty. It can be very 
useful in special cases, but here we want to study the 
question of robustness from another side. How effectively 
can we protect information from noise, against which we 
are unprepared? We can seek answers for this question 
in the following way. We have a channel and a correction 
operation, then we alter the channel a bit, and check 
how effective the unaltered correction operation remains 
on the new channel. In this way we can get information 
on the robustness of a correction operation against com- 
pletely unexpected and subsequently undetermined type 
of uncertainties. 

More formally, let TZq be an error correction operation 
originally used to correct the channel £q. Our aim is to 
characterize the robustness of TZq against the alteration of 
the channel. To alter the original channel, starting from 
£q, we move in the space of channels in a direction given 
by another channel £, in other words we mix the first 
channel with the second one. Mathematically, we define 
channel mixing as the convex combination of channels. 
The operator elements of the mixed channel £' are the 
following: 



{E' J }={Vl 3 7E,}w{^E fc } 



(3) 



where {E^} and {Efc} are the sets of the operator ele- 
ments of £o and £ , furthermore 7 € [0, 1] is the mixing 
parameter. Notation l+J means here that if {E^} and {Ek} 
contains an element which is the same in these two sets 
then this common element will appear two times in {E'^}. 
It is easy to see, that the obtained operator element set 
still represents a quantum channel, as channels form a 
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convex set, so the trace preserving constraint ([2]) remains 
valid. 

During channel mixing we observe the change in the 
effectiveness of TZq which is measured by F c h(TZo o £'). 
However, it is worth comparing it with the best possibil- 
ity allowed by the given encoding and noisy channel £', 
instead of just working with the pure value of it. Ana- 
lyzing whether the effectiveness of TZq remains close to 
the best possible value or not while we alter the original 
channel by mixing it with some £ , it turns out against 
which type of noise is the correction operation sensitive, 
and against which type it is robust. 

Thus the robustness of the TZq correction operation for 
a given £q channel and £ alteration can be characterized 
with the difference between two entanglement fidelity val- 
ues. First, we must determine the entanglement fidelity 
for the case when the effect of the mixed channel £ ' is 
corrected by the best possible TZ* correction operation — 
which gives the biggest possible entanglement fidelity — , 
then we must check, compared to this optimal value, how 
much smaller is the entanglement fidelity we get by using 
the original TZq correction operation: 

AF ch (7) = F ch {TZ* ofy- F ch (K o£' y ). (4) 

This difference can be greater than zero because of two 
reasons. Firstly, it is possible that the examined TZq cor- 
rection operation is not even optimal for the £0 channel. 
The F c h(7?-o £0) — F^TZq o £ ) difference, which is in- 
dependent from the mixing parameter 7, characterizes 
optimality, i.e. how effective the TZq correction for the 
£0 channel is. Secondly, it can be greater than zero, if 
TZq is not robust. We want to use the expression "ro- 
bustness" to characterize the feature of the error cor- 
rection operation, which shows how the efficiency of TZq 
changes with the mixing parameter. Thus it turns out 
that the robustness against the alteration of the chan- 
nel can be measured by the inspection of the 7 depen- 
dence of the AF c h function. If AF c h grows fast with 7, 
then the TZq correction is sensitive to the £ alteration, 
or else when TZq remains closely as effective as it was 
for the £q channel, it is robust. The robustness could 
be defined as the derivative ^AF c h, however it is more 
useful to characterize the robustness with a domain in 
which the efficiency of TZq is not worse than it is for the 
original channel plus a given value. The boundary of the 
(^-robustness domain can be defined as the channel £' 
where 75 is the maximal value, for which it is true that 
AF ch (j)-[F ch (TZ* o£ )-F ch (TZoo£ )] < 6, if 7 < 75. The 
whole (5-robustness domain of a correction operation TZq 
around the channel £q consists of all the channels which 
are inside the boundary. 

To determine the boundary of the robustness domain 
we need to know the second term of |@J, i.e. how good 
entanglement fidelity we get for the channel £' using the 
TZq correction operation. It can be seen, that we do not 
need optimization to determine this term. Using the op- 
erator elements of the mixed channel £!y as in ([3]): 

F ch {TZ o£') = (l- 1 )F ch (TZoo£ )+ 1 F ch (TZoo£) . (5) 



That is, this term varies linearly in function of 7 as the 
distance from the channel £q grows. This follows from 
our definition of channel mixing and the linearity of the 
trace. 

More difficult it is to handle the first term of eq. ([!]). 
For a given channel, the correction operation which gives 
the greatest entanglement fidelity in general can only be 
determined by numerical methods, which should most 
obviously produced as the solution of a semidefinite pro- 
gramming problem. In the literature authors mostly fol- 
low this approach 01, 0, Q , although there are also alter- 
nate attempts 

HE!. Eventually, in the case of a special 
class of channels — the Pauli channels — the optimal cor- 
rection operation can be generated analytically d, [l(| • 
This is a very essential result, as many important chan- 
nels belong to this class, furthermore it also can help 
understand more complex cases, where only numerical 
methods are available. 



III. OPTIMAL QEC FOR PAULI CHANNELS 
A. The standard QEC procedure 

The Pauli group on n qubits is denoted by G n . It is 
generated by all 3n one qubit Pauli operators Xi, Yi, 
Zi acting on the ith qubit. It has order 2 2n+2 , because 
for each i we can have /, Xj, or Zi, plus a possible 
overall factor of ±1 or ±i. If a channel has a set of opera- 
tor elements which contains scaled elements of the Pauli 
group, i.e. {E.; = ^faigi\ where a» > and gi 6 G n , then 
we call it Pauli channel. Because of the trace preserving 
condition of the channel, aj — 1 has to hold. 

The QEC procedure begins with an encoding. This 
can be seen as a mapping Uc : Ti — > K. where TL is 
the Hilbert space of the system which is the source of 
the quantum information — it is spanned by the original 
qubits — , K, is the Hilbert space of a larger system, and 
Uc is the encoding isometry which defines the code space 
C = Uc(W) C JC spanned by the logical qubits. The 
Knill-Laflamme theorem Q states that the channel £ 
is correctable perfectly on C iff there exists a set of op- 
erator elements {Ej} of £ which satisfies the equation 
PcEjEjPc = CiSijPc, where Pc is the orthogonal projec- 
tion onto C and the Cj are scalars. It means that all the 
syndrome subspaces Si := Ej(C) are mutually orthogonal 
and the effect of E, after the encoder can be written as 
EiVc — y/CiAiWiUc where W, is an isometry between 
the code and the syndrome subspaces and Aj is a unitary 
operator. In this case the errors can be corrected by a 
projective syndrome measurement followed by a unitary 
operation which depends on the result of the syndrome 
measurement. Actually this is (A^WiUc)^, the inverse of 
the encoding operation followed by the error operator. 

The efficiency of the QEC depends on the code con- 
struction. Certainly the error operators of a general chan- 
nel do not meet the requirements of the Knill-Laflamme 
conditions. The task is to find an encoding operation 



4 



for which the most probable errors of a given channel 
prove to be correctable. But if the noise affects each qubit 
independently, i.e. the interaction operator between the 
qubits and their environment is tensor product of single- 
qubit operators, furthermore the noise is weak, then the 
most probable errors are those which affect only one 
qubit, therefore the aim is to find codes on which these 
errors can be corrected. We call this strategy in this pa- 
per the standard strategy of QEC. The success of the 
standard strategy does not depend on the actual form of 
the noise. It was the method of the early results of QEC, 
like e.g. the five qubit code [13, HH, and up to now this 
is the main stream of the QEC research. 



B. The optimal QEC construction 

As it can be seen from the previous subsection, after 
fixing the errors required to be corrected and construct- 
ing the code, the correction operator is determined. This 
is quite clear in case of stabilizer codes. The generators 
of the stabilizer group determine not only the code sub- 
space, but also the syndrome subspaces as the intersec- 
tion of the ±1 eigenspaces of each generator and the Aj 
unitaries are determined by the implicit assumption that 
the most probable errors required to be corrected are the 
one qubit errors. 

The case of Pauli errors and stabilizer codes is very 
special because every Pauli group error rotates the code 
subspace onto one of the mutually orthogonal syndrome 
subspaces [l7|. This is because all the elements of the 
Pauli group are either commute or anticommute with 
each other. In the case of an [[n, k]] stabilizer code 
this means that the code subspace C which is stabi- 
lized by the subgroup S — (g%, . . . ,g n -k) < G n is ro- 
tated by the operator g £ G n onto the subspace stabi- 
lized by ((-l)« fll ,...,(-l)*»-* 5n _ fc ) ) where ft = if 
[ffiSi] = and qi = 1 if {g, g{\ = 0. It is convenience 
to label these syndrome subspaces as S q by the indices 
q = q n _ k 2 n - k - i +---+q 1 2°. 

It is clear that the class of Pauli operators which map 
C onto itself form the N(S) normalizer of S which is also 
a subgroup in G n . Furthermore, all the Pauli operators 
from a g-N(S) left coset maps C onto the same syndrome 
subspace S q . Let us fix a representative element from 
each of these cosets and denote them by W q . After this, 
an arbitrary element g of the Pauli group can be writ- 
ten as g = nW q where n £ N(S) and W g are uniquely 
determined. The stabilizer is a subgroup in its own nor- 
malizer, S < N(S). Therefore by fixing a representative 
element i r A p from each of the left cosets of S which are 
contained by the normalizer, i.e. i r A p S C N(S), we can 
write n = i r A p s, where r, p and s £ S are uniquely de- 
termined for each n. Here we write down the scalar fac- 
tor i r explicitly, (r = 0, ...3), which can be done because 
if A P S c N(S) then i r A p S C N(S) too. Actually, the 
A p normalizer operators are the logical Pauli operators 
which act on the logical qubits as the Pauli operators on 



the physical qubits therefore they are self-adjoint opera- 
tors. We conclude that an arbitrary Pauli error operator 
can be written as 

E = = i r A p W q s (6) 

where s £ S, p = 0, . . . , 2 2k - 1 and q = 0, . . . , 2"- fe - 1 
are uniquely determined by E. The index r is unimpor- 
tant, without loss of generality we can assume that it is 
zero because the i r E operators transform the density op- 
erator the same as E, as can be seen from eq. ([T]). So, for 
a given encoder, we classified the Pauli error operators. 
All these are in one of the left cosets of the stabilizer, 
ApW q S, where q means that the error maps C onto S q 
and p specifies a further unitary operation. 

After this classification of the error operators let us see 
the error correction. All of the operators from the same 
coset ApW q S transform the code subspace equivalently, 
that is they cause the same error which can be corrected 
by the same error correcting operation. (When there are 
more than one important error operators in the same 
coset which are required to be corrected, then it is the 
case of degenerate QEC.) Therefore it is convenient to 
reduce the error operators of the form (J6]) with a given p 
and q into one operator ^a p . q A p W q , where a p . q = J^i a ^ 
and the summation index runs over all i for which Ej can 
be written as ([6]) with the given p and q. 

Due to the effect of the error operator a pq A p W q the 
code moves to the syndrome subspace S q and there it 
gets transformed further by the A p operator. Further- 
more, this error occurs with probability a p _ q . Notice that 
if this is the only error operator in the channel (i.e. if 
a p.q = 1), then this channel is perfectly correctable with 
the U^WjAp operator, as this is exactly the inverse of 
the encoding operation followed by the error operator. 
If there is another error operator in the channel, which 
is indexed by a different q integer, then this means that 
with some probability, the code can get into some other 
syndrome subspace, which is orthogonal to the previous. 
This is not a problem, as — according to the spirit of the 
Knill-Lafiamme theorem — after the syndrome measure- 
ment, the code is already projected into one of the S q 
subspaces, so we can apply such a correction operation, 
which has two different operator elements, one corrects 
the first error, and the other corrects the second. But if 
more error operators of the channel move the code into 
the same syndrome subspace in different ways, then — as 
stated by the Knill-Lafiamme theorem — all of these can 
not be corrected perfectly at the same time. It can be 
seen from the following that we can choose any of them, 
and — in case of Pauli errors — can also correct it. But only 
one p for each q. Here it concludes, how effective the cor- 
rection will be for a given channel: the correction will 
be the best, if for each q we choose and correct the most 
probable error. In contrast, the standard QEC chooses to 
correct the few-qubit (more often only the single-qubit) 
errors. 

Now we show, that the intuitive strategy above really 
gives the greatest entanglement fidelity. Let the index of 
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the most probable error operator belonging to each S q 
syndrome subspaces be denoted by p*. Now the opera- 
tor element set of the optimal 1Z* correction operation 



can be obtained in {U^WjA p . } 2 



J C vv g rt P*l 9 =0 

correction operation the entanglement fidelity 



form. Using this 



F ch (TZ*o£) = 

= EE 

p,q q> 

= E a ^> 



Tr (UtW^VWWUe^ 



Wj,W, 



(7) 



5„_„rPr. and 



We used that U C U^ = 

the logical Pauli operators are mutually orthogonal in 
terms of Hilbert-Schmidt scalar product, Tr(PcA p A p /) = 
8p.p' dim(C) = <5 Pj p'2 fe . By this formula, the entanglement 
fidelity in the case of Pauli channels is thus the sum of 



the . 



-p q ,g 



probabilities of the errors selected for correction 



in each syndrome subspace. It is evident, that we get the 
maximal entanglement fidelity, if we design the error cor- 
rection operation to correct the error with the greatest 
probability in each syndrome subspace. Strictly speaking, 
we have proved only that among the (2 2fc ) 2 " different 
QEC operations which can be given in the Kraus form 



-k_ 
o 



there does not exist any better than the one with p q = p* . 
However, there is not any kind of QEC operation which 
is better. A complete proof can be found in [lol ]. 



C. Comparison of the standard and optimal QEC 
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FIG. 1: (Color online) The channel fidelity for the phase- 
damping channel using the code in Table J^a) | and (i) the 
optimal error correction operation: solid curve, third order in 
p; (ii) the standard QEC operation: long-dashed line, second 
order in p; using the code in Table H|T)| and (iii) the optimal 
error correction operation: dotted-dashed line, second order 
in p; (iv) the standard QEC operation: short-dashed line, sec- 
ond order in p; finally the case of no error correction (dotted 
line) is also included for comparison. The optimal correction 
is better on the full domain of the noise parameter < p < 1 
not only than the standard curves but also than the no error 
correction curve. 
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We shortly mention the cases in which the correction 
operation based on optimization could give significantly 
better result than the standard strategy. In these cases 
one should study carefully which are the most probable 
errors. When 

(i) the noise does not act independently on each qubit, 
i.e. the interaction operator is not the tensor prod- 
uct of single-qubit operators; 

(ii) the noise is not weak; 

(iii) even though the channel is tensor product and the 
noise is weak, it can still occur that a single qubit 
error has smaller probability than a several qubit 
error. 

A good example on this third case is the phase damping 
channel, in which X and Y type errors do not occur at 
all. We use the representation of this channel given by 
the Kraus operators 

E =(l-f)/, El .(f)z. (9) 



TABLE I: Generators and logical Pauli operators for the five 
qubit codes. 

Therefore if we use against this noise the five qubit code 
15] which can be defined by the generators from Table 
Hta)| (see [H|), and according to the standard strategy 
we use fifteen syndrome subspaces to correct all the fif- 
teen different single-qubit errors then we get significantly 
worse channel fidelity than in the optimal case, when we 
use the 2x5 syndrome subspaces originally designed to 
correct the single-qubit X and Y errors to correct the 
two-qubit Z errors. Actually, all the two-qubit Z errors 
can be corrected in this way by this five qubit code. The 
result can be seen in Figure [1] where the solid line shows 
the channel fidelity computed by the optimized QEC op- 
eration 1Z* , this curve is cubic in the noise parameter 7 in 
contrast to the long-dashed curve which is only quadratic, 
showing the channel fidelity computed by the standard 
QEC operation. 

Moreover, there is difference between the standard and 
optimal QEC from another point of view. We can also 
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use the five qubit code developed by Laflamme et aLp3 
which can be defined by the generators of Table [ flb)| to 
correct the errors of the p hase dam ping channel. The 
two codes defined in Table [ pa)| and |(b)| are essentially 
the same in the sense of the standard strategy i.e. both 
are [[5,1,3]] codes which correct perfectly all one-qubit 
errors. However, the optimal QEC operations based on 
these codes are very different. As can be seen from the 
dotted-dashed line in Figure [TJ using the code from Ta- 
ble [ yb)| the correction of all two-qubit errors arc not 
possible, therefore the curve is only quadratic. In fact 
the difference of the two codes also can be seen even 
from the curves of standard QEC (short- and long-dashed 
lines), because the two codes correct different subset of 
the higher order (three- and more-qubit) errors. 

Contrary to the above mentioned three cases, the stan- 
dard strategy is acceptable in most cases, when the chan- 
nel is not very asymmetric, i.e. the probability of each 
single-qubit error is about the same. In this case it is easy 
to see that the most likely errors will be the singlc-qubit 
errors, and the traditional QEC correction operation is 
just designed to correct these, thus it is clearly opti- 
mal. A good example on this is the depolarizing channel, 
for which the optimal entanglement fidelity just equals 
with the value obtained from the conventional QEC er- 
ror correction. This result — considered surprising by the 
authors of [6(] — is obvious in the light of the above. 

IV. ROBUSTNESS OF ERROR CORRECTIONS 
ON PAULI CHANNELS 

Now we draw conclusions about the robustness from 
the results of the previous section. For simplicity we as- 
sume that £ , £ and so £' are Pauli channels, so we can 
directly apply the analytical results of the previous sec- 
tion. Therefore the operator elements from eq. ([3]) can be 
written equivalently as 




where a p>q and a Pi9 are the numerical coefficients of the 
reduced operator elements of the channels £ and £, re- 
spectively, and a' p q := (1 — ^f)a p , q + ia>p,q are the nu- 
merical coefficients of the reduced operator elements of 
the channel £'. As can be seen from the formula 
to obtain information about the robustness we need the 
optimal correction for the mixed channel £' . Thus we 
apply the expression ([7]). For this we need to determine 
the most probable one among the errors which move the 
code subspace C into a given syndrome subspace S q , i.e. 
the greatest a Pj(J for a given q. 

To see what exactly happens in each syndrome sub- 
space with the optimal correction operation during mix- 
ing, let us first take a simple case, a Piq = 8 PiP 6 qt q where 
p and q are fixed. In this case only in the syndrome sub- 
space Sq will there be any change in the ratios of the a PA 
probabilities, so there is no need to examine any other 



syndrome subspace, as in those the most probable error 
remains unchanged during the mixing. With the change 
of the 7 parameter in the mixed channel only the prob- 
ability a'p q grows, the other probabilities are decreasing. 
There are three cases: 

(i) cip^q is smaller than some other a P *-,q, i.e. in the 
channel £o it is not the ApW q error operator be- 
longing to the channel £ which is corrected by the 
current optimal correction, TZq. In this case, how- 
ever the probability of the error ApWq is increasing 
as (1 — l)ap,q + 7 while we mix the channels, but 
the optimal correction remains unchanged until the 
decreasing probability (1 — 7)a P |,g remains greater. 

(ii) ap^q is uniquely the greatest error probability in 
the syndrome subspace Sq. In this case it will in- 
crease further during the mixing while the other 
error probabilities will decrease, so ApWq will be 
the most probable error and the optimal correction 
will be unchanged during the whole mixing. 

(iii) Lastly it can occur that a Pj q is one of the great- 
est error probabilities of £q in the syndrome sub- 
space Sq, i.e. there is more than one error opera- 
tor with the same probability. (In other words, at 
least two different indices p* q ^ p" ^ . . . exist, 
for which a p *_.q = a p v> >q — and p^ — p.) In 
this case the optimal correction of the channel £q 
is not unique. We can choose to correct any of the 
most probable errors and we get the same fidelity 
by these completely different correction operations. 
But when the channel is altered by £, for an arbi- 
trary small value of 7 the optimal correction oper- 
ation will uniquely be the one which corrects the 
error A p Wq. 

We can summarize the three cases in the following 
statement. Supposing that TZq — which is designed to cor- 
rect the error A p * Wq — is an optimal correction operation 
for the channel £q. Then it is also optimal correction for 
the mixed channel £' if p ^ p* q and the inequality 

(l-7)%>(l-7)% + 7 (11) 

holds, or if p = p~ and 7 is arbitrary. 

To examine the robustness domain of an error correc- 
tion operation TZq we need to study the fidelities in eq. 
((4]). First we assume that TZq = TZq is the optimal error 
correction of the original channel £q- It does not depend 
on 7 under the conditions stated in the previous para- 
graph, therefore AF c h = as can be seen from eq. ((H). 
Thus, the boundary of the zero-robustness domain is that 
£' channel, for which 70 is the greatest possible solution 
of inequality (fTTj) , 

a p%,q ~ a P,q 
7o = — r , 

or 70 = 1 if p = p*q . For the channel £ ' even the correc- 
tion operation TZq is optimal but not uniquely. Further- 
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more, exceeding this point the optimal correction oper- 
ation changes abruptly and a completely different error 
correction will become optimal. Before this 70 boundary 
point the optimal fidelity goes as 

F ch (Ti; o E'J = 5^(1 - 7 K;,, = (1 - j)F ch (TZ* o fib) , 

1 

which is the same as the linear curve in ([5]), so the 
A-F c h(7) defined by (|3} is identically zero. (We used that 
FchCR-o fi) = 0.) But when exceeds a' * 5 , i.e. in- 
stead of the correction of the error A p ;W^, the correction 
of ApWg will be the optimal solution, and the optimal 
fidelity goes as 

F ch (Ti; o fi;) = ]T(1 - 7 )a p .,, + (1 - 7)0*5 + 7 • (12) 

As can be seen comparing this equation and eq. ([5]), the 
boundary of the (^-robustness domain for a 8 > is that 

fi' channel, for which 

10 

8 



In contrast, for every other TZq correction operation of 
type ([5]) which is not optimal for the fio channel, 75 ~ 8 
without constant term. 

In general when we mix the starting Pauli channel fi 
with an arbitrary other Pauli channel fi which has more 
operator elements belonging to different p, q index pairs, 
then in a given S q syndrome subspace all error probabil- 
ities can vary differently, and this happens in more than 
one syndrome subspace at the same time. In this case 
the originally optimal error correction operation which is 
designed to correct the errors A p *W g in each syndrome 
subspace S q remains optimal also for the mixed channel 
£' as long as the inequalities 

(1 - 7)a p ., g + ia p ^ q > (1 - l)a p . q + 7 a Pi9 (13) 

hold for all p and q. The greatest value of 7 for which 
they are all true defines a 70 for which £' lo is the bound- 
ary of the zero-robustness domain of TZq . Going through 
this point the optimal error correction operation changes 
abruptly and the fidelity goes linearly with 7. Therefore, 
the boundary of the ^-robustness domain for a 8 > is 
that £' l5 channel, for which 75 = 7 o + const-<5. In contrast, 
for every other TZq correction operation of type ijHJ) which 
is not optimal for the fio channel, 75 ~ 8. It is one of our 
main conclusions that the optimal error correction does 
not only give better entanglement fidelity but it remains 
robust in general, except for the case when the greatest 
common solution of eqs. (|13p 70 is zero. 

Using the following geometric picture we can eas- 
ily determine and visualize the boundaries of the zero- 
robustness regions, in which the optimal TZ* correction 
operation is the same. Since all the Pauli channels which 
has the same effect on the code subspace can be given 



by the numbers < a p<q < 1, and because of the trace 
preserving condition J2 P q a p,q = 1; a U °f these channels 
can be represented by a point of a 2 n+k — 1 dimensional 
simplex which has 2 n+k extremal points. By definition, 
we assign the perfectly correctable channels as the ex- 
tremal point at the origin of the simplex. These are the 
channels which can be reduced to A0W0 since Sq = C 
and so Wo = / and also Ao = /. The other extremal 
points of the simplex correspond to the channel reduced 
toA p W g) (p,g)^(0,0). 

A channel fio has unique optimal correction operation 
TZq iff the p* indices are uniquely determined for all q 
by the strict inequalities a p *, q > a p . q for all q and p 7^ 

p*. Mixing fio with another channel fi the optimal error 
correction operation will remain TZq until eqs. (| 13[) hold, 
i.e. a' p , q is maximal among a' p q . The point given by the 
a' p q numbers is on the boundary between two (or more) 
domains if there exist at least one index q and at least 
two indices p* 7^ p** for which a' p , q — a' p ,, q > a' 
where p can be arbitrary. These points are on the 2 n+k — 2 
dimensional plane which is determined by the following 
2 n+k — 1 points: all the extremal points except A p *W q 
and Ap*»W g and the midpoint between these two latter 
extremal points. 



V. DISCUSSION OF PAULI AND NON-PAULI 
CASES 

In connection with Pauli channels and error correc- 
tion operations based on stabilizer codes our main con- 
clusions are the following. In general, the optimal QEC 
operation is robust. It means that the QEC operation 
which is optimal for a channel fio remains also optimal 
in completely unchanged form for other channels fi^ that 
can be given by small alteration of fio. More precisely, 
the channels which has the same optimal QEC opera- 
tion form a zero-robustness domain in the space of all 
Pauli channels. There are finite number of such domains 
and they cover the whole space of Pauli channels. If a 
channel is an interior point of such a domain then the 
optimal QEC operation is robust against any Pauli type 
alteration of the channel, i.e. mixing the channel with 
any other Pauli channel fi, the boundary point fi 7o of 
the zero-robustness domain is in non-zero 70 distance. 
Moreover, the boundary point of an arbitrary 8 > 8- 
robustness domain is outside the zero-robustness domain 
therefore 75 > 70 and lim^o 7<5 > 0. Summarizing, if 
we have a more or less reliable channel model to describe 
the interactions between the system and its environment, 
then choosing the optimal QEC operation we can protect 
against the uncertainty of the channel model without any 
knowledge about it. Only when we are on the boundary 
between two zero-robustness regions is the case worse. In 
this case there are more (at least two) completely dif- 
ferent optimal QEC operations giving the same maximal 
channel fidelity. By choosing any of them we can find such 
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alteration of the channel — actually that one which could 
be corrected by the other optimal QEC operation — for 
which -fs starts to increase linearly with 8. In other words 
there is no such QEC operation which is resistant against 
any type of channel alteration. If we do not know some- 
thing about the uncertainty of the channel, then we can 
not protect against this uncertainty. However, this case 
is exceptional, most of the channels are interior points of 
a zero-robustness domain. 

In contrast, for an arbitrary channel, in the case of er- 
ror correction with any other correction operation which 
has the form ([5]) except the optimal error correction, the 
js value of the ^-robustness region starts to increase lin- 
early with 6. This is also true for any other type of QEC 
operation, having a growth with higher powers of 5 at 
most in special cases. Based on this we finally can con- 
clude that for Pauli channels the optimal 1Z* correction 
operation is better not only in entanglement fidelity but 
also in robustness than any other correction operation, 
assuming the channel is not on the boundary surface of a 
zero-robustness domain. Therefore we should be careful 
when the standard QEC differs from the optimal. The 
latter one takes slight assumptions about the channel, 
but it does not guarantee the robustness at all. Surpris- 
ingly the channel specific optimal solution can give better 
guarantee for the robustness. 

These results are valid for Pauli channels and stabi- 
lizer codes, which is a toy model for the general case, 
but based on some numerical results we believe that the 
main conclusions are similar also in the case of arbitrary 
channels. To demonstrate the general situation we study 
in detail a few examples. For this purpose we choose two 
non-Pauli channel models, the amplitude-damping chan- 
nel, and the less familiar "pure states rotation" channel 
defined by Fletcher et al. in [9(. The amplitude damp- 
ing channel is probably the most well known non-Pauli 
channel. The single qubit operator elements of this chan- 
nel are the following: 



En 







1 



K V^~P, 



E 1= ^ 
\0 j 



This channel model is physically well motivated. It de- 
scribes the decay of a state |1) into the state |0). The 
states |1) and |0) can be excited and ground states of an 
atom, or a one-photon state and zero-photon state of a 
single optical mode, so p is the probability that the ex- 
cited state has decayed to the ground state, or a photon 
has scattered out from the optical mode. 

The pure states rotation channel is in contrast a physi- 
cally much less motivated channel model. It serves for us 
only as an example for a non-Pauli channel besides the 
amplitude-damping, broadening the scope of our results. 
The single qubit operator elements are the following: 



En 



o 

n(V) 



and 



E± =0 



cos (V) sin (f) ±cos(^)cos(f) X 
±sin(^)sin(f) sin (V) cos (f; 



Here a and (3 are determined by the trace preserving 
condition. In this channel the parameter (j) characterizes 
the strength of the noise. Regarding the value of the other 
parameter, for the simplicity, in this paper we use only 
the = f§ value. 

As we mentioned previously, in case of non-Pauli chan- 
nels we have no analytical method to determine the op- 
timal error correction operation. Therefore hereafter we 
follow the literature 0,0, [|[ and based on the observation 
that for a fix encoding it is a special convex optimiza- 
tion problem (a semidefinite program) , we use numerical 
methods instead. In our examples, as encoding operation 
we use the standard five qubit code defined by the gen- 



erators of Table 1(a) 



As a fist step, we study the robustness of the optimal 
QEC operation of a Pauli channel against a non-Pauli 
altering. We have studied four different cases by mixing 
two types of Pauli channels with two types of non-Pauli 
ones. As Pauli channels we use the phase damping chan- 
nel (see eq. (j9])) and the depolarization channel. We use 
the representation of the depolarization channel given by 
the Kraus operators 



En — 



x _3p 
4 



I , E fc = ( - ) a k , k = x,y, 



In the above mentioned four cases we studied system- 
atically the dependence of our results on the noise pa- 
rameters of the two initial channels that were mixed. We 
found that the picture is similar in each case. There are 
no qualitative differences in the results while we modify 
the noise parameters of the depolarizing, phase damping 
and amplitude damping channels between zero and 0.7 
and the parameter <j> of the pure states rotation channel 
between and |j. In Figure [2] we present typical results 
for all of the four cases. In the figures, the first and the 
second terms of the expression AF c h(7) can be seen, that 
is the optimal fidelity ^^(72.* oE ' ) (solid line) and the lin- 
ear term F c h(7\!.n ° £L) (dotted-dashed line), respectively. 
It can be seen that the two curves are almost identical 
at a neighborhood of the initial Pauli channel. (Really, 
close to the Pauli channel, there is no difference in the 
two curves in the order of numerical error.) Neverthe- 
less we can not determine exactly the border of the zero- 
robustness domain, and in principle Sq could be zero, but 
6j starts to increase very slowly with 7, and by no means 
linearly. This region could be the generalization of the 
zero-robustness domain of the pure Pauli case. Increasing 
the distance from the original Pauli channel we leave this 
domain and reach another domain where <5 7 clearly has 
a linear term as a function of 7. Between these domains 
there is a transitional region which could be the gener- 
alization of the boundary surfaces of the zero-robustness 



9 




FIG. 2: (Color online) The two terms of the AF robustness in function of 7 for the mixing of Pauli and non-Pauli channels 
7 = denotes the initial Pauli channel on each figure. The mixtures are 

damping channel (p = 0.3); (b) the depolarizing channel (p = 0.3) and the pure states rotation channel 

0.3) ~ 



the phase damping channel Jp 



the depolarizing channel (p = 0.3) and the amplitude 



5tt 
12 ' 1 



57T 1 

36 ' 



(d) the phase damping channel (p - 



0.3) and the amplitude damping channel (p - 
and the pure states rotation channel (9 = f§, </> = §§ )• The solid line denotes the optimal channel fidelity term, the dashed and 
dotted-dashed lines denote the linear terms for the two initial channels, and the dotted line denotes the approximate derivative 
which is normalized to its maximum. 



domains. To a first approximation we could summarize 
our observations based on the analysis of these fidelity 
curves, that the longer the period which the F^TZ^oS^) 
curve slicks to the linear term, the more we can consider 
the corresponding correction operation robust against the 
given alteration. General robustness domains would be 
those domains, in which the optimal correction operation 
changes slowly with the alteration of the channel, and the 
borders between the regions would be those parts, where 
it changes faster. The necessary condition for the possi- 
bility of robust error correction would be that the channel 
is inside a robustness region. 

To analyze in more detail the above observations and 
to clarify the differences and similarities with the pure 
Pauli case we study the change of the optimal QEC op- 
eration with the mixing parameter 7. In the Pauli case 
it was completely unchanged in the zero-robustness do- 
main but it was changed abruptly at the boundary of 



the domain. In the general case, to measure how fast is 
the change in the optimal QEC operation TZ* we com- 



pute the || ^X; I 



derivative of the Choi matrix XX which 



is the matrix representative of the superoperator describ- 
ing the QEC operation TZ*. Moreover, || A|| := ^Tr(AtA) 
denotes the Hilbert-Schmidt norm. We compute this 
derivative approximately, it is plotted by dotted line in 
Figures ffia)fl(d")j We observe that TZ* varies very slowly 
near the Pauli channel, i.e. if 7 is small. Even more in- 
terestingly we can see that there are abrupt changes in 
7?.* , denoted by the peaks of the dotted curves. The bet- 
ter the numerical precision, the higher these peaks are, 
indicating discontinuous jumps in the TZ* optimal QEC 
operation. Between these non-continuous changes TZ* is 
not constant but changes relatively slowly as a smooth 
function of 7. This indicates that the optimal QEC op- 
eration changes in two different ways. The first case is 
similar to the domain boundary of the Pauli case, where 
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FIG. 3: (Color online) The two terms of AF in function of 7 
for the mixing of the amplitude damping channel (p = 0.5) 
and the pure states rotation channel (8 = §f ,0 = §J), and 
the approximate derivative j|^X^||. Notations are similar to 
those of Figure [2] 



FIG. 4: (Color online) The enlarged version on Figure 2(c) 
in the small vicinity of 7 = 1. The second term of AF for 
the non-Pauli side was subtracted from the first, thus the 
red line shows the AF function. The dotted line is again the 
normalized approximate derivative ||^X*||. 



in a syndrome subspace one of the error operators be- 
came more likely than the others. We believe that this 
similarity can not be accidental, however at this moment 
we can not see how it could be precisely generalized to 
the non-Pauli case. 

We also have to mention that there are more break- 
points in the optimal fidelity curve which can not be seen 
because of the scaling of the figure. However, changing 
the scale and enlarging the vicinity of a point where TZ* 
is discontinuous, we always observe a small breakpoint 
in the optimal fidelity curve. It is important to note that 
between these breakpoints, where TZ* varies analytically, 
the optimal QEC operation is robust in the sense that 
small change in 7 causes a change of only second or- 
der in AF. But in the non-analytic points this is not 
true, the optimal TZ* is not uniquely defined and there is 
no good choice against an arbitrary channel alteration. 
These points however form a discrete set in the 7 € [0, 1] 
interval, and a subset of zero measure in the space of all 
possible channels. In the light of the above results we can 
assert that the optimal QEC operation of a Pauli channel 
is also robust against non-Pauli channel perturbations. 

On the other hand, the curves of Figure [2] can also 
be viewed from the other side, i.e. from the 7 = 1 end- 
point. This point describes the initial amplitude damp- 
ing or pure states rotation channels, respectively. Going 
backwards from this point in 7 we can study the robust- 
ness of the QEC operation which is optimal for a non- 
Pauli channel against a Pauli type alteration. As can be 
seen from ([5]), the second term of eq. (j4]), F^TZq o £'") 
varies linearly in function of 7 for an arbitrary channel 
and an arbitrary alteration, i.e. as the distance form the 
channel £0 grows. This is indicated by the dashed line 
in Figure [2] In contrast to the case of Pauli channels we 
can see from the dotted lines that the derivative of the 
Choi matrix is not zero. This means that the TZ* optimal 
QEC operation changes by arbitrary small alteration of 



the channel. As a consequence Sq — 0, i.e. there is no 
finite zero-robustness domain around a non-Pauli chan- 
nel. However, from the fidelity curves we can see that 
the optimal term F^TZ* o^) slicks to the linear term 
F c h(TZo o £') and 75 starts to increase only with higher 
powers of 8 and not linearly. Therefore the optimal QEC 
operation proved to be more robust than any other QEC 
operation, for which the 75 value of the (5-robustness re- 
gion starts to increase linearly with <5. Only when we 
are in a point where the optimal QEC operation changes 
abruptly — at the peaks of the dotted line — is the case 
worse, very similar to the boundaries of zero-robustness 
domains in the case of Pauli channels. Wc can sec an in- 



teresting example in Figure 2(c) Very close to the ampli- 



tude damping channel there is a very narrow peak in the 
derivative indicating an abrupt change in the 

optimal QEC operation, therefore the optimal and the 
linear lines split at that point. Enlarging the vicinity of 
the amplitude damping channel (see Figure^]) we observe 
that the optimal and the linear terms are almost iden- 
tical when 0.9982 ;$ 7 and the difference between these 
terms grows linearly as 7 decreases, i.e. the distance from 
the initial amplitude damping channel is increasing. The 
picture is very similar to the Pauli case, nevertheless the 
robustness domain is very small. 

To make sure that our statements are general, we also 
studied a case where two non-Pauli channels are mixed. 
For this purpose we mix the amplitude damping and 
the pure states rotation channel. Again, we scan the 
p £ [0, 0.7] and <f) £ [0, §§] noise parameter intervals of 
the channels. As a typical example we plot the p = 0.5 
and 4> = || case in Figured Our observations are essen- 
tially the same as in the previously discussed case. Ac- 
tually, we can also see such non-Pauli-non-Pauli chan- 
nel mixing cases anywhere in the figures of Pauli non- 
Pauli cases, when we take two channels corresponding to 
nonzero 7 points and mix them. 
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The results stated above strongly indicate — however 
not prove — the following consequences. In contrast to the 
pure Pauli case, zero-robustness domains, where the op- 
timal recovery operation is not changes at all, does not 
cover the whole space. Instead, TV can change in two dif- 
ferent ways. The first is similar to the case of boundaries 
of Pauli zero-robustness domains. In these special points 
71* is not unique, it changes abruptly and the fidelity is 



non-analytic function of the mixing parameter. Between 
these points TZ* changes analytically with the mixing pa- 
rameter like the fidelity, and therefore the optimal QEC 
operation can said to be robust. 

A more precise study of the distribution of the non- 
analytic points in the space of all quantum channels is 
the task for the future. 
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